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BODIES DURING SHCfRT-TIME, HIGH-SPEED FLIGHT 
By Hbu Lo 


SUMMARY 


A short and simple method is presented for the determination of 
transient skin temperature of conical bodies for short— time, high-speed 
flight. A differential equation is established for this purpose, 
giving the fundamental relations between the transient skin temperature 
and the flight history. For the heat-transfer coefficient and 
boundary— layer temperature, which are needed in the differential 
equation, Eber*s experimental results for conical bodies under super- 
sonic conditions are adapted and summarized in a convenient way. The 
method is applied first for flight at constant altitude to illustrate 
the effect of acceleration on transient skin temperature. It Is then 
applied to arbitrary flight. Several examples are given; for one 
example measured data are available and are in good agreement with the 
calculations. 


INTRODUCTION 


When air flows over a body the air immediately adjacent to the 
body is brought to rest by skin friction. As a result the air is 
heated to a higher temperature and hence there is heat exchange between 
the air and the skin. This phenomenon is generally termed as 
"aerodynamic heating." 

At high speed the temperature Increase of the air is very large 
and the aerodynamic heating problem becomes of great concern to 
designers. The problem is related to the characteristics of the 
boundary layer and the local heat— transfer coefficient. In reference 2 
Wood gave a method for the determination of the skin temperature at 
supersonic speed, using formulas for heat— transfer coefficient and 
boundary-layer temperature, derived for flat plates at subsonic speed. 
In reference 2 Scherrer made a more theoretical approach for the 
determination of skin temperature of a body of revolution at supersonic 
flight. Both papers, however, dealt with equilibrium skin temperature 
for steady flight conditions at constant altitude. 
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In several German papers (reference 3, for instance) it was shown 
that, for a short-time flight during which the speed and altitude vary 
with tine, the transient shin temperature may he considerably lower 
than the equilibrium skin temperature. In the present report, 
therefore, emphasis is given to the transient skin temperature, rather 
thaw the equilibrium skin temperature. A differential equation is 
presented for this purpose, and for the heat-transfer coefficient and 
boundary-layer temperature needed in the differential equation, Eber*s 
experimental results (reference 4) for conical bodies under supersonic 
conditions are adapted and summarized in a convenient form for 
immediate application. If, however, better experimental data become 
available, they can be adapted readily to the present method. 

To show the effect of acceleration on skin temperature the 
simpler problem of flight at constant altitude is treated first. 

More general flights are then discussed with several examples. In 
one example, data obtained from the Naval Research Laboratory, 
Washington, D.C., giving skin temperatures for V2 missile number 21 
fired on March 7, 19^7 at White Sands, N. Mex., in connection with 
upper atmosphere research (reference 5), are used for comparison with 
the calculated results. 
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SYMBOLS 

acceleration, ft/sec^ 
recovery factor 

specific heat of air at constant pressure, Btu/lb/°F 
specific heat of skin material, Btu/lb/°F 

function of T sk and implicitly of t (see equation ( 16 ) ) 

acceleration due to gravity, taken as 32.2 ft/sec^ 

heatr-transfer coefficient, Btu/(sec)(sq ft)(°F) 

reference heat— transfer coefficient corresponding 

to B = jt/6 (or 30°), Z = 1 ft and H = sea level, 
Btu/(sec)(sq ft)(°F) 

thermal conductivity of air, Btu/( sec )( ft)(°F) 
characteristic length, ft 
time, sec 
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time when certain velocity is. reached by uniformly 
accelerated body starting from rest, sec 

correction factor for nose angle (see equation ( Ih- ) ) 

correction factor for characteristic length ( see 
equation (lk)) 

correction factor for altitude (see equation ( lk) ) 

heat absorption capacity of skin, Btu/sq ft/°F 
altitude, ft 

mechanical equivalent of heat, taken as 778 ft— Ib/fctu 
heat flowing into skin due to skin friction, Btu/sec/sq ft 

heat lost through radiation, Btu/sec/sq ft 
Reynolds number 

free— stream temperature of air, °F abs. 
skin temperature at time t a , °F abs. 
equilibrium skin temperature, °F abs. 
initial temperature, °F abs. 
boundary-layer temperature, °F abs. 
skin temperature, °F abs. 
stagnation temperature, °F abs. 
velocity of flight, ft/sec 

specific weight of skin material, lb/cu ft 

total apex angle of conical nose, radians 

emissivity, assumed to be O.k in numerical, examples 
of this report (e = 1 for perfect black body) 


* 
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k 

0 1 free-etream density of air, lb/ou ft 

p free— stream density of air at sea level, taken as 

bL 0.07657 lb/cu ft 

11 coefficient of viscosity of air, lb/sec/ft 

T Bkin thickness, ft 

ANALYSIS 

Fundamental Equations 


The skin tsmperature of a body can be determined from a considera- 
tion of the amount of heat flowing into the skin and' that lost by the 
skin. The balance of the two is the heat absorbed or given up by the 
skin from which the change in skin temperature can be calculated. 

Heat flowing into skin .- For heat flowing into the skin during a 
short-time, high-speed flight, the most important factor to be con- 
sidered is the "aerodynamic heating." Experimental results indicate 
that the heat flowing into the skin due to skin friction Q^_ can be 

determined from the following empirical formula: 

= h^T BL - T ak ^ Btu/sec/sq ft (l) 

where the boundary— layer temperature T BB and the heat— transfer 

coefficient h have been studied experimentally by Eber (reference k) 
and will be discussed later. 

Heat received by the skin from other sources, like solar 
radiation, radiation from surrounding atmosphere, heat from interior 
or other parts of the body, and so forth, are not considered in thiB 
report. 

Heat los t by the skin .— The heat lost oy the skin can be (1) heat 
loss through radiation, (2) heat loss through artificial cooling, 

(3) heat loss to interior or other parts of the body. The heat loss 
through radiation Qg can be determined from the Ste fan-80 ltzmanri 

formula: 

$2 = ! *.8 x 10 ~ ^-3 eTg^ Btu/sec/sq ft (2) 


where the value of the emissivity e depends on the surface condition 
of the 3kin. For a perfect black body, 6=1. 
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The heat loss other than radiation to the exterior will not be 
considered in this report. 

Heat balance equation .— Equation (1) gives the rate of heat flow 
into the skin and equation (2) gives the rate of heat loss by the skin. 
The difference of the two (Qi — is the heat left to heat the 

skin, Btu/sq ft/sec. During an interval dt, the total heat to be 
absorbed by the skin is therefore 

(Ql — Q2) < ii' Btu/sq ft 


If the temperature rise of the skin during this interval dt is dTgj^, 
the heat absorbed by the skin can also be expressed by 


G dTg^ Btu/sq ft 


where G, the heat absorption capacity of the skin, Btu/sq ft/°F, 
is the product of the specific heat of the skin material c 8 , the 

specific weight of the skin material W, and the skin thickness t . 
From data in references 6 and 7, representative values of G for 
steel and alloys of aluminum and magnesium at 520° F are obtained 
and plotted against skin thickness in figure 1. At higher temperature 
the values of G are increased at the rate of 0.256 t * per 100° F 
for steel and 0.062Vr* per 100° F for aluminum where t* is the 
skin thickness in inches. For magnesium alloy the temperature effect 
can be neglected. 

Equating the preceding two expressions for the heat absorbed by 
the skin and substituting and Q2 from equations (1) and (2) 

gives 


dT sk 

G “dT + hT 


a k 


+ 1^.8 x 10 13 eT 


sk = hT BL 


(3) 


which is the basic equation for transient skin temperature. 

The simplified equation .— Equation (3) is a nonlinear differential 
equation with variable coefficients. A simplified equation can be 
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obtained if the radiation loss can he neglected. The simplified 
equation is 

cLT 

«-# - “sk * “%L C*) 


As will he shown later (example 3) the radiation loss ordinarily 
contributes only a small portion to the transient skin temperature 
for short-time, high-speed flight. For such flights the use of 
equation (4) will not introduce any appreciable error, while a 
great deal of time and labor can be saved in the computation work. 

Equilibrium temperature .— If a body flies with a constant speed 
at a constant altitude for a sufficiently long time, the skin is 
heated to a temperature such that the heat flowing into the skin 
per second is for practical purposes the same as the heat lost, and 
the skin temperature becomes constant. The skin temperature then is 
called the equilibrium temperature T e , corresponding to that speed 

and altitude. The equilibrium temperature T e can be determined 

dT sk 

from equation (3) by dropping out the first term since ^ = 0. 

Thus, 

4.8 X 1O“ 13 cT 0 4 + hT e = hT BL (5) 


If the radiation loss can be neglected, the equilibrium skin tempera- 
ture will be the same as the boundary— layer temperature . 

In equations (3)* (4), and ( 5 ) the two parameters h and Trt. 

are needed before the equations can be solved. In the following, 
Eber*s experimental data on these two parameters will be discussed 
and summarized. 


Eber*s Experimental Results 

In 1941 Eber (reference 4) made a series of wind-tunnel tests on 
cones of various vertex angles to determine the boundary— layer 
temperature and heat-transfer coefficient h at high speeds. 

The results he obtained are adapted in this report for the determina- 
tion of skin temperature for conical bodies in flight. These results 
will be discussed in the following and will be summarized in a 
convenient form with certain simplifications . 
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Boundary— layer temperature, Tg^.— When the air stream is brought 

to rest isentropically, the temperature of the air is called stagnation 
temperature. When the air is brought to rest by skin friction, the 
process is usually not isentropic and the temperature of the air is 
lower than the stagnation temperature . In the latter case, if there 
is no heat flow across the skin, the temperature of the air 
immediately adjacent to the skin is called the. boundary— layer 
temperature Tg L and is found to be closely related to the stagnation 

temperature. In the actual case where there is heat exchange between 
the air and the skin, the boundary— layer temperature Tgg is only an 

artificial term used in the empirical equation (l) and may not 
necessarily be realized at any point in the actual boundary layer. 

The stagnation temperature Tg.p can be calculated theoretically 
from the following statement of Bernoulli's equation: 



( 6 ) 


where and T]_ are the velocity and the temperature of the free 

air stream. Values of T]_ at various altitudes are given in 
reference 8 and are replotted In figure 2 of this paper. 

For constant value of Cp, equation (6) becomes 


t ST 


- Ti = x 


IV 


2 Jgc_ 


( 7 ) 


For the actual case, however, Cp varies with temperature. It has 

been pointed out by Wood (reference 1) that at high velocity a 
correction for variable c p will lower the stagnation temperature 

considerably (about 20 percent lower at Mach number 8). An exact 
solution of Bernoulli's equation for variable Cp will give a set 

of stagnation temperature curves for various altitudes and Mach 
numbers. In this report, however, it is chosen to plot the stagnation 
temperature rise (t st - T^ against free-stream velocity Vj_, 

instead of the stagnation temperature against Mach number. The 
result is a single curve for all altitudes. This curve, for which 
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the derivation is given in appendix A, is plotted in figure 3, together 
with the stagnation temperature rise for Cp = 0.24 Btu/lb/°F. For 

low speeds, no correction for variable c p is necessary. 

The boundary-layer temperature • T BL is related to the stagnation 

temperature by an empirical factor c, called recovery factor. The 
recovery factor c is defined as 


c 


Tbl 


t st ” T 1 


( 8 ) 


Eber*8 experimental results show that c varies with the total 
vertex angle 3 of the cone but is practically independent of 
velocity. The variation of c with 3 is not very large. For 3 
ranging from 20° to 50°, an average value c =0.89 can be used 
with a maximum error of about 2 percent. Therefore, 

( t BL ~ T l) « 0.89 (t st - T x ) ' (9) 

Since the stagnation temperature rise (t&t — can be 

taken as a function of only (fig. 3)» the boundary— layer-temperature 

rise can also be taken as a function of V^j that is, a single curve 

for all altitudes. This curve is shown in figure 4. For 

Cp = 0.24 Btu/lb/°F, the boundary— layer-temperature rise can be given 

by the following, expression: 



( 10 ) 


Heat— transfer coefficient, h .— The heat— transfer coefficient h 
can be determined from Eber*s empirical formula which is taken from 
reference 4: 


h = (0.0071 + 0.0154 \/3^)|(R) 0 * 8 


where the Reynolds number R, as defined by Eber, is equal 
to (piV.IAi). Equation (ll) can be rewritten as 


( 11 ) 
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h = (0.0071 + 0.015 b\TB) 



( 12 ) 


In the derivation of equation (ll), Eber used the boundary— layer 
temperature as the reference temperature for values of k and p, 
the thermal conductivity and coefficient of viscosity of air, 
respectively. For the characteristic length l, Eber used the 
entire length of the surface of the cone. Although there are 
different opinions regarding what is the correct length to be used 
as characteristic length, when Eber*s experimental results are applied 
to actual flight conditions (for instance, in reference 3 Kraus and 
Hermann used half the total length of the cone surface as charac- 
teristic length) Eber's original definition will be followed in this 
paper. The value of h thus obtained represents an average value. 
Instead of local value, of the heat— transfer coefficient. 

A simplified method for the evaluation of the factor (k/p 0 * 8 ) is 
now given. For a given velocity and altitude, the boundary-layer 
temperature can be determined from figures 2 and k, and values of k 
and p corresponding to this temperature can be obtained from 
reference 9 for temperatures below 21+00° F absolute. This vas done for 
various velocities at three different altitudes ( sea level, 100,000 ft 
and 190,000 ft) and the ratio k/p®* 8 was calculated and plotted 
against velocity in figure 5* For boundary— layer temperature higher 
than 2k00° F absolute the curve is shown dotted and is obtained 
by extrapolation. 

In figure 5 the two curves corresponding to 100,000 feet 
and 190,000 feet form two limiting values of k/p®. For any other 
altitudes from sea level to 370,000 feet, approximately, values 
of k/p^* 8 against velocity will fall within these two limiting 
curves. Since both limiting curves do not differ very much from 
the sea— level curve, it is Justified to use the sea— level curve for 
all altitudes up to 370,000 feet. Therefore, k/p®* 8 is a function 
of velocity only. 

In connection with this simplification, it should be kept in 
mind that values of k and p at high temperature are obtained by 
extrapolation and any effects of change of air composition at high 
altitudes on k and p are not considered. 

Equation (12) is now reduced to four factors which are functions 
of 3, l, altitude and velocity, respectively. For convenience of 
computation, a reference heat— transfer coefficient h', corresponding 

to 3 = A (or 30°)» 1=1 foot, and at sea level, is confuted from 

equation (12) and the results are plotted in figure 6. This reference 
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heat-transfer coefficient h* is a function of velocity only. For 
other nose angles, characteristic lengths and altitudes, the heat- 
transfer coefficient h is simply the reference heat— transfer 
coefficient h* multiplied by the three correction factors Fg , Fj, 

and Fg . Thus, 


where 


h = F p FjF H h» 


(13) 


0.0071 + 0.01^4\/T 
0.0071 + 0.0154\T«^> 


F Z " j0.2 



(14) 


Values of Fg, Fj, and Fg are given in figures 7* 8, and 9. Values 
of (pi/Psl) > where pg L is the free-etream air density at sea level, 

were taken from the NACA standard atmosphere table ( reference 10) for 
altitudes below 65,000 feet and from tables V(a) and V(b) of 
reference 8 for altitudes above 65,000 foot. 

Application of Eber*s results. — In the foregoing, Eber*s 
experimental results on and h are represented by simple curves 

as functions of velocity and altitude. For any prescribed flight path 
where the velocity and altitude are given as functions of time, values 
of Tjjl and h can be expressed as functions of time. Table 1 is 

prepared for this purpose and the operations in table 1 are self- 
explanatory. Knowing and h as functions of ' time, one can 

solve equation (3) or (4) for the transient skin temperature, and 
equation (5) for the equilibrium skin temperature. 

Certain facts should be borne in mind, however, in the application 
of Eber*s work, particularly to high-altitude flight conditions. First 
Eber*s experimental results were obtained over a limited range of 

Reynolds number, 2 x 10^ to 2 x 10^j for flight conditions where the 
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actual Reynolds number falls outside of this range (which is usually 
the case for flight at high altitudes), extrapolation is necessary. 
Second, Eber*s experiments were carried out at low altitudes. At 
high altitudes. where the air density is very low, the flight may 
enter into the slip— flow domain where the recovery factor, and' 
possibly the heat— transfer coefficient, may be greatly affected. 

(For illustrations, calculations were made to determine the Mach 
numbers and Reynolds numbers from the flight trajectory of the 
V2 missile used in example 1 of this report. The results of these 
calculations are indicated. by points in figure 10, where the curve 
dividing the two domains is taken from reference 11. For this 
particular example the flight enters the slip— flow domain as soon 
as an altitude of approximately 130,000 feet is reached.) Third, 
in the evaluation of k and p, the effects of change of air 
composition at high altitudes are neglected and values of k and p 
at temperature greater than 2*4-00° F absolute are obtained by 
extrapolation. Finally, the properties of the atmosphere at high 
altitudes are taken from the tentative tables of reference 8. With 
all these uncertainties the application of Eber*s results to high 
altitudes may introduce a large percentage error. However, since 
the heat— transfer coefficient at high altitude is small, the effect 
on skin temperature is not large. 


Solutions of Equations 


Solution of equation ( 3 ) Equation (3) is a nonlinear equation 
with variable coefficients, of the first order but the fourth 
degree. If equation (3) is written in the following form 



(15) 


where 


h_ 4.8 x 10~ 13 e T *4 

G i BL G i sk 5“ i sk 


(16) 


it is readily recognizable that equation (15) can be Bolved con- 
veniently by Runge and Kutta*s numerical method (reference 12) which 
is summarized in appendix B. In table 2 Runge and Kutta*s method is 
arranged in a suitable way for the solution of equation (15). A \ 
convenient time interval At is first chosen. The smaller At is, 
the more accurate the results will be . The operation of this table 
starts on the first line, proceeding from left to right, and then 
the succeeding lines. The function f, corresponding to t Tg^ 

at its left, can :be obtained from equation (l6). The final results 
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of table 2 give the skin temperature at the end of the interval At. 
If the table is repeated, the skin temperature at the end of 2At, 
3At, and so forth, can be obtained. If the time interval chosen 
is not small, a correction can be made as explained in appendix B, 
where an illustrative example is also given. 

Solution of equation (IQ .— Equation (4) is a linear differential 
equation of the first order provided the heat absorption capacity G 
is considered to be independent of skin temperature. The general 
solution is 



where D, the constant of integration, can be determined from the 
initial condition 


T . = T, at t = 0 
sk 1 


The parameters h and Tg^ in equation (17) usually cannot be 

expressed in simple analytical terms. Therefore, equation (17) has 
to be integrated numerically. Table 3 is provided for this purpose. 
The time and labor required for carrying out the computations in 
table 3 are much less than those required for table 2. 

Solution of equation (5) .— Equation ( 5 ) is a simple quartic 
algebraic equation. By use of Eber*s experimental results on h 
and T bl , equation ( 5 ) can be solved for the equilibrium 

temperature T e for any altitude and velocity combination. For 

illustrative purposes, the variation of equilibrium temperature 
with velocity for several altitudes is shown in figure 11. The 
value of 6 used is 0.4. 


RESULTS AMD DISCUSSION 


The method for the calculation of transient skin temperature, 
as discussed in the preceding sections, is applied first to flight 
at constant altitude to illustrate the influence of acceleration on 
skin temperature , and then to arbitrary flight conditions. 
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Flight at Constant Altitude 

For Illustration of the effect of acceleration, the simplified 
equation (4.) will he used to determine the skin temperature for the 
following example . 

Assume a body, starting from rest, is accelerated uniformly 
to 5000 feet per second. It then maintains this speed until the 
equilibrium temperature is reached. Subsequently, the body decelerates 
uniformly to zero velocity again. How does the skin temperature change 
with time during these three periods of flight? The body h as a conical 
nose angle of 3O 0 , a conical surface length of 1 foot, a skin heat- 
absorption capacity G - 0.6, and is traveling at a constant altitude 
of 50,000 feet. 

Period of uniform acceleration .— The velocity of the body at any 
instant during this period is given by the following equation 

= at 

p 

where a = acceleration, feet per second . The Initial condition 
is Igk ;s T^ at t s '0. 

Since the velocity is given as function of time, T-g^ end h 

can be determined from table .1 and the skin temperature from table .3, 
all as functions of time. The results are given in figure 12 fpr 
five different values of acceleration: a = 2g, 5g, lOg, 50g, 

and 00 (for infinite acceleration the temperature— time curve is but 
a single point) . The dotted curve in figure 12 gives the skin tempera- 
ture T a at the time t a when the body reaches 5000 feet per second . 

The larger the acceleration is, the lower the skin temperature will be . 

In figure 13, the skin temperatures for a = .2g and lOg are 
plotted against V^, together with the curve showing the variation 

of equilibrium temperature with velocity,. Since radiation loss is 
neglected in this case, the equilibrium temperature is the same as the 
boundary— layer temperature . The temperature difference between the 
"transient" and the "equilibrium" curves is the " temperature lag," 
which is greater for larger accelerations, as shown in figure 13. 

Period of constant velocity .— If the body maintains its speed 
at 5000 feet per second after it reaches this velocity, h/G 
and TgL are no longer functions of time, and the solution of 

equation (17) becomes 
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T ek - TBL + De - ^ 0 ^ for t S t a (18) 

where D, the constant of integration, is determined from the 
condition 

. T sk - T a when 1 m *a : ' . 

Here T a , the skin temperature at time when the "body first reached 

the velocity of 5000 feet per second is different for different, 
accelerations . 

Equation (18) is solved. for the five different accelerations 
a = 2g, 5g> 10g, 50g, and «. Results are plotted in figure 14. 

(The temperature curves before the body reaches 5000 ft/sec are taken 
from fig. 12.) The skin temperature approaches exponentially the 
equilibrium temperature T e , in this case equal to Tg^ corresponding 

to 5000 feet per. second, at the altitude .of ,50, OCX) feet. 

Period of uniform deceleration .— If after the body maintains 
5000 feet per second for a long time, it starts to decelerate, the 
velocity at any instant is given by 

= 5000 + at 

where t starts when -the body begins to decelerate. The initial 
condition is 

T sk = T BL at > = 0 

Table 1 and table 3 again can be used to solve equation ( 17) and the 
results are plotted against velocity in figure 15 for two different 
decelerations a = -2g, — lOg. Again the equilibrium temperature 
curve is' also given. The transient skin temperature in the case of 
deceleration is higher than, the equilibrium temperature throughout 
the velocity range. The temperature lag is therefore on the adverse 
side. • . • •' 


General Flight Conditions 

In the general case, the body is changing its altitude as well 
as its velocity. Three examples are given in the following. In all 
cases., the simplified equation (4) is used for the calculation of 
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transient skin temperatures. In example (3) the more accurate method 
is also used and the relative importance of radiation Iobb is 
discussed . 

Example 1 .— The V2 missile number 21 fired on March 7, I9L7 at 
White Sands, N. Mex., had a conical vertex angle of 2 6°, a conical 
surface length of approximately 7 feet and a skin of 0.109— inch steel 
after a certain distance from the nose.. The above specifications 
and the flight path of V2 as shown in figure 16 were obtained from the 
Naval Research Laboratory, Washington, D.C. The boundary— layer 
temperature and the heat— transfer coefficient are determined by 
means of table 1 and are plotted against time in figure 17. The 
skin absorption capacity G is obtained from figure 1 to be O.U76 
corresponding to 520° F absolute, which is conservative. The initial 
skin temperature is known to be 5U6 0 F absolute. Table 3 is used 
to determine the skin temperatures and the results are plotted 
against time in figure 18. 

In figure 18 the measured skin temperature for the steel skin 
of the same missile is also shown. The measured data are obtained 
from the Naval Research Laboratory. The agreement between the 
calculated and measured results is good. 

Example 2 .— An arbitrary velocity and altitude diagram is 
assumed, as shown in figure 19, including descending path as well as 
ascending. The missile is assumed to have a nose angle of 30°, a 
characteristic length of 4 feet, and a skin heat-absorption capacity 
of 0.6 Btu/sq ft/°F. The boundary— layer temperature T BL and heat- 

transfer coefficient h are determined by use of table 1 and are 
plotted against time in figure 20. Table 3 is then used to calculate 
the skin temperatures. The results are plotted in figire 21. The 
skin temperature during descent becomes higher and higher because of 
the greater density at the lower altitudes and the higher velocity 
as it comes down. 

Example 3 «— To determine the relative importance of the 
radiation loss, the calculation of skin temperature for the missile 
in example 2 for the first 80 seconds is repeated except now the 
more accurate method is used where the radiation loss is not 
neglected (the emissivity e is assumed to be 0.L). Table 2 is 
used for this purpose and the results are plotted in figure 22 
together with the results obtained from example 2 where the radiation 
loss is neglected. The discrepancy between the two is very small, 
approximately 3° F. 

In fact, for most missiles the radiation loss plays only a Bmall 
part in the determination of skin temperature and the simplified 
method can be used to great advantage in saving time and labor. 
Unfortunately, there is no simple criterion to predict when the 
radiation loss can be neglected. Generally speaking, if the skin 
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heat capacity is not too small, say not below 0.5 Btu/sq ft/°F, the 
8 kin temperature is not expected to be higher than 1000° F absolute 
and if the radiation loss is not the dominating factor for a long 
period of time, the radiation loss can be neglected. If the final 
skin temperature is completely unknown, it is advantageous to start 
with the simplified method. From the maximum skin temperature obtained, 
an estimate of the radiation loss and its effect on the skin temperature 
can be quickly made. For instance, the maximum skin temperature of the 
missile of example 2 is about 615° F absolute. Conservatively assume 
this is the skin temperature for the entire period of 80 seconds. 

The change of skin temperature due to radiation loss during this 
period is then 


at* - s ( 4 - 8 x 10-136 1 T sic U ) 4 

=* oTS^* 8 x 10-13 x °*^ x 6l 5 U ) x 80 


= 3-7° P 


which is negligible. A similar calculation for example 1 indicates 
that for the first 65 seconds of flight the radiation loss affects 
the maximum skin temperature approximately 6° F. 

For flight conditions where the radiation loss is the dominating 
factor for a long period, as from 80th second to 220th second in 
example 2, the radiation loss should be investigated. During this 
period, h is zero and equation (3) beooues 



4.8 



= 0 


The general solution is 



(1 


14.4 x 10' 


,- 13 , 


t 

G 


+ 



( 19 ) 


where D, the constant of integration, can be determined from the 
condition 
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L sk 




80 


when t = 80 


where 


( T 4o is the skin temperature at the 80th second. At end 


of 220th second the skin temperature can he calculated from 
equation ( 19 ) and is found to he lowered only 6° F. 


CONCLUDING REMARKS 


1. A differential equation taking into account aerodynamic 
heating and hody radiation is established for the calculation of 
transient skin temperature for any prescribed flight history. Runge 
and Kutta*s numerical method is recommended for the solution of the 
differential equation. 

2. A simplified differential equation which neglects body 
radiation is also given, and can be used in many cases to great 
advantage in saving both time and labor. 

3 . Eber*s experimental results on the boundary— layer temperature 
and heat-transfer coefficient, to be used in the differential equation 
are summarized in a convenient way for immediate application. Tables 
and charts are also provided to facilitate the solution of the 
differential equation. 

4- . The calculated skin temperature for a V2 missile Is In good 
agreement with the measured data. . 

5. The heat-absorption capacity of the skin has an important 
influence on transient skin temperature. The heat-absorption capacity 
is greater and consequently the temperature lag is larger if (a) the 
skin is thicker, (b) the material is denser, or (c) the specific heat 
Is higher. 

6 . When the air is heating the skin the temperature lag due to 
the heat capacity of the skin is in the favorable direction; that is, 
tends to lower the skin temperature . When the air is cooling the 
skin, the temperature lag is in the adverse direction; that is, tends 
to keep the skin at high temperature. 

7 . Eber*s experimental work was conducted under certain limited 
conditions (short testing time, small temperature difference, limited 
Reynolds number, and so forth). More refined experimental values for 
a wider range are desirable. 
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8 . Because the atmospheric properties at extremely high altitudes 
as given in reference 8 are tentative, and also Because at high 
altitudes the flight actually enters the "slip-flow" domain, a closer 
investigation of the problem at high altitudes is needed. 

9 . For more accurate results, investigations of heat exchanges 
other t han those considered in this paper are necessary. 


Langley Memorial Aeronautical Laboratory 

Rational Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX A 

STAGNATION TEMPERATURE RISE FOR VARIABLE c p 


In the following, a single curve for all altitudes is obtained 
for the calculation of stagnation temperature rise at various 
velocities Vi when the specific heat c p is considered to he a 
function of temperature. The Bernoulli equation can he written as 


V dV + JgCp dT = 0 


(Al) 


Integrate equation (Al) between state 1, the free -stream condition, 
and state 2, where the velocity is zero. 

0 P t ST 

V dV = / JgCp dT 

1 U T 1 
or 


2 


l ST 

Jg I Cp dT 
T, 


(A2) 



Values of the integral in equation (A2) can he obtained from 
table 1, reference 9 (see also page 58, reference 9) for given values 
of T]^ and Tg.j, and the velocity V^ can be computed. A set of 

curves can thus be obtained. 

However, if it is chosen to plot ^Tgip — T^ against V^ 
with T-l as parameter, the set of curves will practically fall into 
one single curve for all values of T^ ranging from 392° F absolute 

to 63O 0 F absolute, corresponding to the minimum and maximum free— 
stream temperature for altitudes from sea level to about 370,000 feet. 
In figure 23 two curves are shown, representing the two extreme 
conditions. It is Justified, therefore, to use a single curve for 
the calculation of stagnation temperature rise. 
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The boundary— lay e r— tempe ra ture rise ~ can obtained 

by multiplying the stagnation temperature rise by the recovery 
factor c. The result will be a single curve similar to the curve of 
stagnation temperature rise but with ali the ordinates decreased by 
the ratio c. Figure b shows the curve of boundary— layer— temperature 
rise at various speeds V^, corresponding to c = O. 89 . 
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APPENDIX B 

’ t - • f I , 

RUNGE AND EOTTA'S METHOD OF NUMERICAL INTEGRATION 


Runge and Kutta’s method of numerical integration can be applied 
to differential equations of the following type/ provided that the 
initial. condition ie known. ' 

j| = f(x,y) (Bl) 


The initial condition is 

4 * l • « » ' ■ *, ! * 

y =? y Q when jt • x*, , (B2) 


The derivation of Runge and Kutta'B method can be found in 
reference 12. The following table is provided for the calculation. 


• X 

y 

f ' 

fx(Ax) 

Operation 

x o 

fo 

f (xoiyo) 

*1 


X 4^ 

x o*T 

V§»i 


^2 

q 2 +q 3 » 

V* 

36^2 

f^b+^,y 0 +|q^ 

^3 

Sum » • 




% 

q=^ Sum = 


X 1 = x© + Ar Xi = y Q + q 


For convenience of application to the present problem, Runge 
and Kutta*s method is arranged in the form of table 2 of this paper. 

A suitable interval At should be chosen. Starting from the pair 
of initial values t 0 and T Q , one can obtain the values t® and T* 

at the end of the interval At, by carrying out the operations in 
table 2. The operations proceed from left to right of the first line 
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and then the succeeding lines. The function f(t,T 8lc ), corresponding 

to t and T to its left on the same line, can be obtained from 
equation (l6). If the pair of values t' and T* are used as 
initial values and the table repeated, another pair -of values t" 
and T" can be obtained, corresponding to the skin temperature at 
the end of 2At. Table 2 can be repeatedly used in this manner 
until the desired values of T are obtained. 

For illustration, table 2 will be used to calculate the skin 
temperature of the missile in example 2. Assume the skin temperature 
at the end of 20th second is known to be 519° F absolute, and the 
8 kin temperature at end of 24 th second is to be determined. 

Choose At = 2 seconds . 


t 

T sk 

h 

t bl 

G 

f 

fx At • Operation 

20 

519 

0.0103 

544 

0.6 

O.403 

q x =0.806 |(<li+<l4.) 

» 

1.618 

21 

519.4 

.0106 

568 

.6 

.834 

<12=1 *668 ( l2 + ^3 


3.324 

21 

519.8 

.0106 

568 

.6 

.828 

q^=l .656 Sum 

= 

4.942 

22 

520.7 

.0107 

590 

.6 

1.215 

4k “2 .430 q 

= , 

1.644 

22 

520.6 

.0107 

590 

.6 

1.217 

q x =2 .434 ^ (q.i + ik) 

= 

3.299 

23 

521.8 

.0108 

619 

.6 

1.728- 

q 2 *3 .456 q 2 +< i3 

= 

6.886 

23 

522.3 

.0108 

6i9 

.6 

1.715 

q^=3 .430 Sum 

s 

IO.185 

24 

524.0 

-.0109 

640 

.6 

2 .082 

q^ “4 . l64 q 

= 

3.395 


2b 524.0 


Therefore, the skin temperature at the end of 24th second is 
524° F absolute. Notice that a constant value of G is used in the 
above computation. For a more accurate analysis the value of G 
should be based, on T s1e given. 

To improve the result, the following method,- as given by Runge 
and Kutta, can be used for corrections. Instead of At = 2 seconds, 
table 3 is repeated with At = 4 seconds. 
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t Tgjj. h T bl G f f x At Operation 

20 519 0,0103 544 0.6 0.403 q. 1 =1.6l2 1^+^) ■ .4. 976 

22 519.8 .0107 590 .6 1.230 <12=4.920 q 2 +q 3 = 9.704 

24 521.6 *0107 590 .6 1.196 <1^=4.784 Sum = 14.680 

24 523.8 .0109 640 .6 2.085 q^=8.340 q = 4.893 

24 523.9 ; " 

The correction is then given by ' 

1 =. 75(524 - 523.9) = 9.007 

•• • • • •' 15 •• ;,V 

and the corrected skin temperature at the end of 24 th second 1 b 

* - ‘ « * 1 . 1 ‘ t ’ 

Tg k « 524;+ | = 524.007 

Runge and Eutta's methodcan also be applied to simultaneous 
equations of the type (Bl) or to differential equations of higher 
order. For details, see reference 12 . 
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Table, / 

Variation of Tgi_ and h with time 
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t 

sec 

H,rt 

Given 

y: £L 

1 sec 
Given 
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7% 7 

Fh 

fyio 
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1 

l 

i 

! 


' 



^5^ 

p 


Given jh- total nose angle = 

l = nose surface length t ft. - 
fg ( figure 8 ) = 

Fj, ( figure 9) = 











Table 2 

Numerical •solution crF more accurate transient shin 

temperature equation 
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If t\ 7' are used as 'initial conditions and the above table 
repeated ' the skin temperature corresponding to t = t a i 2 At 
can he obtained ' 
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Table 3 

Numerical solution of the simplified 
transient skin temperature equation 


n 

n 

/ 

l 

3 

4 

e/fc. 

m 

t 






ii 

h (Table. 1 ) 






© 

h/G ' 
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msmmm 

0 






e® 






© 

T bl ( Table 1) 






m 







©■ 

WHMSMStflnBKik 

— 






©*@ 
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m 

“5® + Ti 

Ti 






Tsk = ©/© 







6 = heat absorption capacity of skin, ^ 
T = initial temperature of skim, °F abs 








0 


0.04 


o. zo 


0.08 


0.1 Z ' 0.16 


Skin thickness, in. 


Fiqure /. - Heat absorption capacity 
of various metal skins at room 
temperature ( 520 °F abs ) . 
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Free- stream temperature, T t} °F abs 

Fiqure 2.- Free -stream temperature at 
various altitudes ( data obtained 
from reference 3 ) . 
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Free-stream velocity , [// , ft /sec 

FtyureZ.- Stagnation -temperature rise for 
both constant and variable Cp , 



Boundary-loyer-temperature rise, 
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0 2000 4000 6000 8000 10,000 

Free-stream velocity , V t , Ft /sec 

Figure 4. - The boundary-layer- temperature 
rise at various velocities (variable ) w 









, Btu 

Reference heat-transfer coefficient, h ,(sec)($c)ff)(‘F) 
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0 2000 4000 6000 8000 / 0,000 
Freest ream velocity , V t , ft /sec 


Figure 6, - l/onotion of reference heab- 
tro refer coefficient with velocity . 
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Toto / vertex angle, fZ, deg 


Figure 7 . - Correction -factor for fZ . 



Characteristic length, l, ft 
Figure 8. ~ Correction factor for l . 
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Altitude , H , -ft 

(a) From sealevel to / 60,000 ft , 
Figure 3, - Cor reef ion Factor For attitude . 
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Figure 10 r Different domains of flow conditions entered by V2 
missile 2i , fined on March 7> 194 7, 
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4,000 

3,000 

1,000 

1,000 

0 

Figure //. ~ Fouijibrium temperature at 
J various velocities and altitudes * 
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$ gpj Q ( y s ± ( 9Jnj.ojdduuaj. uiyg 







Figure / 2 . - Time history of skin temperature -for uni 
accelerated body at constant altitude , 



Skin temperature, T s p. y °F obs. 
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Free-streom velocity, v, } ft/sec 


Figure !3.~ Comparison of transient and equilibrium 
temperature for uniformly accelerated 
body at altitude 50,000 ft . 
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^qoj 0 ' y s j_ ‘ sunjcusduiaj. u/ys 



^9° J Q ‘ ‘ 9jn+DJ2diU9q_ u/o/q 
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N 



*0 



/(jure 14.- T/me history of uniformly accelerated 

body at alt dude 50,000 ft. Constant Fiyure / 5.- Sk'm temperature of uniform 
Velocity after the body reaches 5000 decelerated body at constant 

ft /sec, altitude Sop 00 ft . 
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tj sec 

Figure / 8 . - Comparison of calculated and measured 
skin temperatures for VZ ItZI , fired on March 
7, 19 4-7 ( measured data obtained from 
Naval Research Laboratory ). 
















Skin temperature } Tsk, °Fabs. 


46 NACA RM No. L7K17a 



.... t , sec 

Figure 22, - Comparison of simplified 
method and the more 
accurate method in 
example 3 . 




]/, y ft /sec 

Figure 23.- Stag nation -temperature rise 
variable c n 









